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A multiscale modeling strategy is developed to study grain boundary fracture in 
polycrystalline aluminum. Atomistic simulation is used to model fundamental nanoscale 
deformation and fracture mechanisms and to develop a constitutive relationship for 
separation along a grain boundary interface. The nanoscale constitutive relationship is then 
parameterized within a cohesive zone model to represent variations in grain boundary 
properties. These variations arise from the presence of vacancies, intersticies, and other 
defects in addition to deviations in grain boundary angle from the baseline configuration 
considered in the molecular dynamics simulation. The parameterized cohesive zone models 
are then used to model grain boundaries within finite element analyses of aluminum 
polycrystals. 


I. Introduction 

P HYSICS-based modeling of fracture begins at nanometer dimensional scales in which atomistic simulation is 
used to predict the formation, propagation, and interaction of fundamental damage mechanisms. These 
mechanisms include dislocation formation and interaction, interstitial void formation, and atomic motion such as 
Coble creep. The development of these damage mechanisms progress into microscale processes such as local 
plasticity and small crack formation. Ultimately, damage progression leads to macroscopic failure modes such as 
plastic yielding of components and large cracks exhibiting Mode I, II, and III opening behavior. Thus, multiscale 
modeling strategies are needed for the purpose of developing a unified description of the hierarchy of processes that 
govern fracture. 

One strategy for multiscale modeling is to map an atomistic molecular-dynamics (MD) simulation of physical 
behavior at the nanoscale into continuum finite element representations for microscale analyses. These damage 
mechanisms can include the physical processes that govern intergranular and transgranular damage initiation and 
propagation in metallic micro structures and may be used to define effective traction-displacement relationships for 
modeling fracture at larger length scales. Once the traction-displacement relationships from MD are known, they can 
be used to define effective cohesive zone elements that can be placed along all possible fracture interfaces in a finite 
element mesh. Finite element simulations can then be used to study the failure properties of the material at larger 
scales. The methodology is very generic and can be used for any interface that is subjected to fracture provided that 
the atomistic simulation captures all relevant mechanisms. However, because of limitations on the domain size in 
the molecular-dynamics simulation (with typical dimensions on the order of tens of nanometers), this direct 
approach may be limited to simulating micro structures with characteristic dimensions on the order of microns. 
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Several examples include fracture along grain boundaries in metals (intergranular fracture), fracture along weak slip 
planes (transgranular fracture), and fracture through interfaces between grains and small second phase particles. 

The key to this approach is the recasting of the atomistic results within the framework of a continuum cohesive 
zone model. Cohesive zone models (CZMs) approximate traction-displacement relationships along a material 
interface and are frequently used in conjunction with the finite element method (FEM). They are embedded into 
zero-thickness CZM elements, which may be placed between continuum finite elements to simulate interface failure. 
In fracture studies at the micromechanics level, CZM elements can be used to predict transgranular fracture if they 
are placed adaptively between continuum finite elements within grains or to predict intergranular fracture if they are 
placed along grain boundaries (GBs) as shown in Figure 1. 

Parameterizing the relative influence of various atomistic mechanisms and crystallographic orientations is 
important if CZMs are to be used effectively in a multiscale modeling strategy. In the present study, molecular 
dynamics (MD) simulations based on a bicrystal model are used to extract the combined elastic and plastic 
properties of a common grain boundary in aluminum under Mode I opening. The information is then recast in the 
form of bilinear CZMs to represent the mechanisms involved in fracture. This basic functional form is then used 
parametrically to study the effect of CZM parameters (representing variations in properties of individual grain 
boundaries) on grain boundary fracture within aluminum poly crystals. 

In this paper, multiscale modeling is developed to study intergranular fracture in aluminum with cohesive zone 
models implemented to bridge length scales. First, the basic formulation, numerical features and parameterization 
of the CZMs based on MD simulation are discussed. Second, the MD analysis, including the extracted mechanisms 
of atomic- scale deformation and fracture, is reviewed. Third, the method for parameterizing the traction- 
displacement response along the fracturing grain boundary is presented. Finally, sample results based on this 
parameterization for both a bicrystal and a poly crystal are discussed. 



Figure 1 : Embedding CZMs along GBs to study micro structural fracture. 


II. Finite Elements Incorporating Cohesive Zone Models 
A. Cohesive Zone Models 

Numerous CZM formulations have been presented in the literature (Needleman, 1990, Tvergaard and 
Hutchinson, 1992; Costanzo and Allen, 1995; Camacho and Ortiz, 1996; Klein and Gao, 1998; Turon, et al. 2004). 
A general review of various forms of CZM is found in Chandra, et al. 2002. Central to most formulations is the 
requirement that the area under the CZM traction-displacement curve represents the work of separation required to 
open a crack or the fracture toughness and is given by Tvergaard and Hutchinson (1992) as 

G c = Jtc/A (1) 

0 
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where G c is the work of separation, x is a general form of the traction, A is a general form of the displacement, and 
is the critical displacement at which complete separation has occurred and the tractions are zero. 

A mixed-mode CZM developed by Turon et al. (2004) has been used in the current investigation wherein the 
normal and shear components of the traction and displacement are interpolated based on mode mixity. A brief 
outline of the basic formulation given by Turon et al. (2004) is presented here. The interpolation of CZM 
components is depicted in Figure 2. An interpolated CZM assumed in a bilinear form in which relative 
displacements between the upper and lower nodes, A = b top - h bo \ indicate the onset of damage, and final failure of 
the cohesive zone is shown in Figure 3. 




Figure 2. Interpolation of CZM components. Figure 3. Bilinear cohesive zone model. 

(From Turon et al 2004.) 

An equivalent relative displacement is defined as 

X. = a /(A 3 ) 2 +A 2 , w (2) 

where the Macaulay bracket is defined as <x> = 1/2 (jc + |x|), A 3 is the non-negative relative displacement in Mode I, 
and Ashearis the combined relative displacement in Mode II and Mode III given by 

A s w=V(Aj+(A 2 ) 2 (3) 


The mode mixity ratio is defined as 


p- 


A 


shear 


^ shear (^ 3 ) 


( 4 ) 


The effect of combining individual fracture modes on fracture evolution is characterized by assumed interactions of 
the various fracture toughness measures, G. The strain energy release rate, G, is computed by relating the critical 
fracture toughness measures for Mode I, II, and III given by G Ic , G IIc , and G IIIc , respectively, together with the 
current measures for the strain energy release in shear, G shear, and the total strain energy release rate, G tot . An 
effective measure for G is given by the model of Benzeggagh and Kenane (1996). This model is utilized in the 
present study and is given by 


G = G k: + ( G IIc - G h 


G \ 

^ shear 
Gtot j 


( 5 ) 
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where G = Gj+ G s h ea r and G shear = G n + G m . The ratio of the combined shear and the total strain energy release rate 
may be expressed in terms of the mode mixity ratio as 


g G shear P 

G m 1 + 2(3 2 - 2(3 

The mixed-mode damage onset displacement is calculated as 

and the mixed mode displacement corresponding to complete failure of the cohesive zone is given by 

A/ =7^ hc+(Gnc-G Ic M_ 

K A 

Assuming a bilinear CZM as shown in Figure 3, a damage parameter, d , is computed as 



The resulting constitutive law is given by 
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( 6 ) 


(7) 


( 8 ) 


(9) 


( 10 ) 


For the bilinear CZM presented in Turon et al. (2004), k is the initial slope of the traction - displacement law, which 
in this case when the GB width is non-negligible, represents the GB stiffness. 

B. Parameterizing Cohesive Zone Models 

While CZMs are a useful method to simulate material failure, because the parameters that are typically used to 
define them are taken from bulk material properties (Iesulauro, 2002), they do not describe adequately the actual 
physics of micro structural fracture. Macroscale values of strength and toughness that are typically input to the CZM 
represent the aggregate responses of thousands or millions of grains, grain boundaries, and defects within the 
specimens from which they were obtained. Thus, these macroscale values do not represent the unique response of a 
particular interface at which a local fracture event might occur. If the microscale predictions are to become 
quantitative, consideration of the local nanoscale properties is required. One possible means of making this 
connection is to use the results of MD models as input to the CZM. This connection would allow more realistic 
simulations leading to accurate predictions of the failure properties of a large class of materials and microstructures, 
even when experimental data is not available. 

Attempts to extract relevant parameters for the decohesion law of a CZM from atomistic (molecular-dynamics, 
molecular-static, or ‘first-principle’) simulations have been made by various groups in the last few years (Gall et al., 
2000; Komanduri et al., 2001; Spearot et al., 2004; Raynolds et al., 1996). The approach in all of these works is 
based on simulating the debonding of a flat interface under a constant tensile strain rate perpendicular to the 
interface. In these references, the system size is between 4 and 8 nm, and the dynamics of the atoms is severely 
constrained by the boundary conditions, which do not allow for Poisson lateral contraction and shear deformation. 
As a result, plastic processes, such as dislocation slip, are strongly suppressed. Consequently, the simulated 
mechanism for interface decohesion in these references reproduces the process of atomic adhesion (strength) rather 
than that of fracture at the interface. 
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The boundary conditions at which a typical CZM element operates in a large-scale finite element (FE) model are 
very different from the ones used in the referenced MD and first principles simulations. Typically, the CZM 
elements are embedded within a system of finite elements that reproduce the elastic and plastic response of the 
surrounding material to both the external load and the crack-tip stress. In contrast, decohesion parameters, such as 
peak stress and opening displacement of the CZM curve, are typically extracted from an atomistic volume of order 
10 nm in each direction in the referenced simulations. The lack of an adequate surrounding volume of material and 
the stiff boundary conditions suppresses the plastic processes, such as dislocation nucleation, limiting the 
accommodation of deformation at the interface and forcing the configuration to debond in an unnatural manner. The 
periodic boundary conditions usually applied in these models cause the simulation to create a response of an array of 
repeating units with a strong overlap of image elastic forces rather than the response of a single specimen unit. 
Consequently, the resulting decohesion curves cannot be directly applied to derive the constitutive laws for CZM 
elements. 

Thus, an MD-based CZM decohesion law for intergranular fracture under local conditions similar to those 
experienced by the CZM element in a polycrystalline FE model is developed in the present work. The CZM 
decohesion law reflects the response of the CZM element to an approaching and propagating crack (Costanzo and 
Allen, 1995; Davila, 2001). The present MD model is a model of crack propagation rather than of adhesion and is 
used to simulate a crack propagating through a flat high-energy grain-boundary in aluminum (Yamakov et al., 
2005). 


III. Development of a Cohesive Zone Model Based on MD Simulation 

A. Configuration Considered in Molecular Dynamics of Grain Boundary Separation 

For metallic micro structures, the effective GB thickness is small. Figure 5 shows a typical 299 GB formed 
between two grains in an fee metal. Common neighbor analysis (CNA) (Honeycutt and Andersen, 1987; Clarke and 
Jonsson, 1993) is used to identify atoms in different crystallographic states: fee (small dots), hep (triangles), and 
non-crystalline atoms (large dots). Atoms with more than 1/3 of their nearest neighbors missing are identified as 
surface atoms (squares), indicating existing vacancies in the GB. The length scale is in units of the lattice constant of 
Al, a Q = 0.405 nm (Mishin et al., 1999). Here, the specific problem of intergranular fracture is considered wherein 
MD simulations of fracture along GBs are used to obtain the atomistic response of crack propagation under applied 
loads along a characteristic length of the GB. 

[7 7 fO] 



Figure 5. Atomistic snapshot giving the crystallography and structure of the GB interface. 

B. Plastic Contribution to the Stress Field Near the Crack 

The results of a molecular dynamics simulation developed by Yamakov and coworkers (Yamakov et al., 2006) 
that reveals many details of the crack propagation and dislocation response in the bicrystal system are summarized 
here. Figure 6 shows the deformation mechanisms and crack growth of the atomic system resulting from four 
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different initial hydrostatic stresses: a = 3.5 (a), 3.75 (b), 4.0 (c), and 4.25 GPa (d). CNA is used to identify atoms 
in different crystallographic states: fee (in gray), hep (in red), non-crystalline atoms (in blue), and surface atoms (in 
green). Thus, a number of different formations are indicated in the figure as follows: © - GB interface; © - twin 
boundary; © - core of a partial or twinning dislocation; © - nanovoid at the crack tip; © - slip dislocation; © - GB 
dislocation; and © - secondary slip. Due to the crystallographic orientation of the grains, the crack induced along 
the grain boundary propagates dissimilarly in the +x and -x directions. The crack tip is blunted in the -x direction as 
is characteristic of ductile fracture whereas the crack tip remains sharp in the +x direction as is characteristic of 
brittle fracture. The deformation in the -x direction is dominated by the creation of a series of nearly symmetrical 
twins. Although twinning is not common in Al, the preference for twinning to dislocation slip is in agreement with 
the Tadmor and Hai criterion (Tadmor and Hai, 2003). 



Figure 6. Snapshots of cracks that have propagated in the MD system. 


C. Development of Traction-Displacement Relationships 

Through the definition of a cohesive zone volume element (CZVE, Yamakov, et al. 2006) - an atomistic analog 
to a continuum CZM element - the results from MD simulations may be recast to obtain averaged continuum 
traction-displacement relationships that represent the cohesive interactions along some characteristic length of 
material. In Figure 7, stress and opening profiles are extracted along the crack growing in the system prestressed at 
4.25 GPa hydrostatic load and exhibiting the deformation shown in Figure 6d. The insert in Figure 7 shows the 
schematic diagram of the slicing of the system volume in the MD simulation and defining the representative regions 
for extracting the parameters for the CZM elements suitable for use in a continuum simulation. 

In Figure 7, a curve fit of f m is shown as a close approximation to the c f^(x) values extracted from the MD for 
the +x crack tip and are consistent with the predominantly brittle fracture seen in this direction. In contrast, the 
plastic deformation processes, including twinning and dislocation emission, that dominate near the -x direction crack 
tip, have a pronounced effect on the stress distribution. The c f^(x) values extracted from the MD for the -x crack tip 
are nearly constant with position, consistent with a strong plastic contribution. This asymmetric behavior is well 
explained by the Rice criterion (Rice, 1992) for cleavage vs. dislocation blunting transition at the crack tip. 
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Figure 7. Stress and opening profiles extracted along the crack growing in a system 
prestressed at 4.25 GPa hydrostatic load. 

These differences in behavior suggest that two separate cohesive laws are needed to describe brittle and ductile 
crack propagation along one GB. The o^(jc) and X(x) profiles in Figure 7 are for a single instant of time of the 
crack propagation. If the entire simulation of the interface debonding is divided into N t equal intervals of time t q 
(q€z[ 0, AJ), many such profiles can be taken of many CZVEs placed along the GB. When plotted in a cfyy vs. X plot, 
each (&' yy (x p ,t q ), X(x p , t q )) couple represents a point cf s y y(k(x p ,t q )). After sorting these data points in order of ascending 
X, so that X( < X i+ i, and taking a moving average (or a consecutive mean): 


A)= 


(2 M + 1) 


M , . 


j—M 


( 11 ) 


in which the results are averaged over M points backward and M points forward from X h a construction of a 
statistically representative traction-displacement function r(X) can be made. The calculated values for x(Z) for the 
brittle and ductile crack tips are shown in Figure 8. Complete details of the modeling, fracture mechanics and 
extraction of the CZM response are given in Yamakov et al. (2006). 



Figure 8. Brittle and ductile traction-displacement relationships from MD simulation. 
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IV. Parameterized Approach to Multiscale Modeling 

Although numerous forms of cohesive zone model (CZM) have been developed in the literature (Chandra et ah, 
2002), a general bilinear CZM is sufficient to approximate the traction-displacement response shown in Figure 8 for 
Mode I crack propagation along the 299 GB configuration. By adjusting the parameters of the bilinear CZM, a 
range of peak tractions and opening displacements can be considered to represent variations in GB decohesion 
properties transitioning from the more brittle to the more ductile type of decohesion. These variations may arise 
from different GB misorientation angles, defining different types of GB structures (Wolf, 1990) as well as the 
presence of vacancies, interstitials, and other defects in the GB interface. The CZM then provides a numerically 
efficient, but qualitative, means of considering a broad range of GBs without atomistic modeling of all of the 
mechanisms and configurations that contribute to the variations in GB properties. 

The traction-displacement relationships, along with corresponding works of separation (Eq. (1)), may take the 
CZM forms presented in Figure 9. The functional form of the CZMs may be varied to represent process zones that 
range from perfectly brittle to perfectly ductile. An intermediate model expressed as a general bilinear traction- 
displacement relationship is also shown in Figure 9 and represents the common situation in which a mixture of 
elastic and inelastic deformation mechanisms exists and characterizes crack propagation. Here, the parameters x° dh 
x° TO , x° br represent the peak tractions of the perfectly brittle, intermediate, and perfectly ductile CZM models, 
respectively. Similarly, A c dt , A c wx , and A c br represent the maximum displacements of the perfectly brittle, 
intermediate, and perfectly ductile CZM models, respectively. 

T 


A 

Opening Displacement 

Figure 9. Possible traction-displacement relationships for perfectly brittle, perfectly 
ductile, and mixed elastic/plastic crack tip behavior. 

V. Microstructural Configurations 

Figure 10 shows the two configurations considered in the finite element analysis implementation of the 
parameterized cohesive zone models. The size of both configurations is chosen to be l x = l y = 380 nm, so that the 
characteristic length of the cohesive zone element placed along the crack path is / x /200 = 1.9 nm, i.e., equal to the 
size of the CZVE used in the MD simulations (Yamakov et al, 2006). In both cases, 6-node triangular finite 
elements are used to represent the elastic continuum having the elastic properties of Aluminum (E=72 GPa, v=0.32). 
Figure 10a shows a simple bi-crystal configuration with a single interface between the two crystals and an initial 
crack of length a Q = 7.6 nm. The model is loaded with a linearly increasing tensile displacement in the y-direction 
that increases from U y = 0 at l x = 0 to U y = 6% strain at l x = 380 nm. Figure 10b shows a more complex configuration 
consisting of a tessellated polycrystalline micro structure having many grain boundary interfaces but with no pre- 
existing crack. The polycrystalline micro structure is simulated under a uniform tensile strain U y loaded 
incrementally to 22% strain. 
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a) Edge cracked tension model b) Random microstructure model 

Figure 10. Cracked tension model for Mode I opening. 


VI. Results for Parameterized CZM 

Parametric studies were performed to assess the effect of the relative amount of elastic and plastic deformation at 
the crack tip on the fracture characteristics of the configurations shown in Figure 10. This was achieved by varying 
the shape of the CZMs using the MD extracted cohesive zone curves shown in Figure 8 as a baseline. Two sets of 
parameterized CZMs were generated as shown in Figure 11. In the first set (Figure 11a), the effect of the relative 
amount of elastic and plastic deformation was varied, resulting in different values of the energy of crack face 
separation or toughness, G, while keeping x mx constant. In the second set (Figure 1 lb), the ratio of elastic to plastic 
deformation was varied by changing the maximum opening displacement, A mx , for CZMs having the same 
toughness. Keeping G constant results in a varied maximum peak traction x mx . In this way, the effect of the shape 
and the area of the CZM curve on the fracture properties of the interface can be demonstrated. In Figure 11, CZMs 
A1 and B1 correspond to the MD-generated results shown in Figure 8. 




a) Series A: Equal peak traction b) Series B: Equal toughness 

Figure 11. Cohesive zone models used in two series of parametric studies. 

The parameters used to define the models shown in Figure 11a and Figure lib are shown in Tables la and lb, 
respectively. 
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Table la. CZM parameters used to simulate different relative elastic and 
plastic behavior at the crack tip. Fracture toughness is a function of the 
increase in plasticity. The peak traction is held constant at ^ mx =5.0 GPa. 


Cohesive Zone Model 

A mx 

G (J/m 2 ) 

Al 

2.17 

5.42 

A2 

2.89 

7.22 

A3 

4.33 

10.84 

A4 

6.50 

16.25 


Table lb. CZM parameters used to simulate different relative elastic and plastic behavior at the 
crack tip. Fracture toughness is held constant at G = 5.42 J/m 2 . 


Cohesive Zone Model 

x°mx (GPa) 

A 0 (nm) 

A c mx (nm) 

B1 

5.00 

1.07 

2.17 

B2 

3.75 

0.80 

2.89 

B3 

2.50 

0.53 

4.33 

B4 

1.75 

0.37 

6.19 


Figures 12 and 13 show the results for the bicrystal configuration of Figure 10a. The length of the crack was 
determined in two ways (see the inset in Figure 12a). First, the length contributions of all CZM element integration 
points that exceeded the peak traction x mx were summed to give the position, a x , of the front edge of the process zone 
at the crack tip where ‘softening’ starts. Second, the length contributions of only those CZMs that exhibited 
complete loss of stiffness were summed to give position, a 2 , where the process zone ends and a free surface is 
formed. In this way, the fracture process zone is determined as the ‘softening’ region between a x and a 2 . 



(b) 

-G = 5.42 
-G = 7.22 

-G = 10.84 
G= 16.25 


a) Measures ai and a 2 given separately 


b) Mean crack length 


Figure 12. Relative crack length vs. percent edge strain for CZM Series A. 


Figure 12a shows the position a x and a 2 of the front and end of the crack process zone as a function of a linearly 
increasing far-field strain. Here, CZM Series A was considered. In all cases, the a x is greater than a 2 as it indicates 
the beginning of the processing zone. Note that the size of the process zone, estimated as the distance between the a x 
and a 2 lines for each G, increases with increasing G. In addition, the two sets of lines exhibit a systematic trend 
towards shorter crack lengths as the fracture toughness, G, increases. The shortening of the effective crack length, 
taken as the mean of the two positions a x and a 2 with increasing G, is observed in Figure 12b. In spite of the fact that 
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the peak traction x mx , corresponding to the onset of damage, was kept constant, the increase of G leads to the crack 
opening at larger strains. 



Percent Edge Strain Percent Edge Strain 


a) Measures ai and a 2 given separately b) Mean crack length 

Figure 13. Relative crack length vs. percent edge strain for CZM Series B. 

Figure 13a presents the analogous results for CZM Series B. The same effect of increasing the process zone is 
found, but this time it is due to the increased size of the CZM plastic region by decreasing x mx while G was kept 
constant. The effective crack length versus applied strain is shown in Figure 13b. Two different regions of crack 
growth can be seen. First, at the moment of opening, the crack is unstable and the mean crack length exhibits an 
instantaneous growth until it reaches a stable equilibrium. The zone of unstable growth is indicated by the shaded 
region in Figure 13b. Second, after reaching the equilibrium state, the equilibrium crack lengths for d\\x mx follow the 
common equilibrium master curve. This result shows that when G is constant, the equilibrium crack length for a 
given strain remains the same and does not depend on the peak traction x mx . Thus, the peak traction in the CZM 
affects only the initial strain of crack opening. 

The results presented in Figure 14 demonstrate that the process zone in both CZM Series, A and B, 
systematically increases when the CZM parametric curve extends towards larger crack opening displacements, 
either through increasing of the fracture toughness at constant peak traction, or through decreasing of the peak 
traction while keeping the fracture toughness constant. Note that the decrease of x mx has a more pronounced effect 
on the process zone than increasing G and indicates that lowering the traction barrier for debonding results in a 
larger increase of the processing zone in front of the crack tip compared to the case when G is increased. 



Figure 14. Process zone vs. l/x mx for Series B and Process zone vs. G for Series A. 
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Although the MD results have only examined Mode I opening under very specific conditions for brittle and 
ductile crack tips along a 299 grain boundary in pure aluminum, they provide a physical, albeit limited to nanometer 
scales, constitutive relation of decohesion that can lead to the construction of a database of CZMs for studying the 
failure of poly crystal micro structures. For example, by using the MD parameterized CZM traction-displacement 
relation in Figure 8a, a sequence of boundary opening states leading to the formation of a dominant microcrack in 
the configuration shown in Figure 10b can be illustrated as shown in Figure 15. 



Figure 15. Process of dominant microcrack formation in a polycrystal model using cohesive zone models 


An examination of the decohesion processes exhibited in the simulation presented in Figure 15 shows the initial 
opening of GBs under normal applied displacements. The CZM for Mode I opening obtained from MD analysis at 
this length scale has an initial GB stiffness that is of similar magnitude to the stiffness of the adjoining grains. Thus, 
these interfaces are shown to open at low levels of applied strain (Figure 15a). At higher applied strains, the 
energetics of GB opening dictate the location at which a dominant crack begins to form (Figure 15b). At a critical 
applied strain, the local coalescence of GB opening forms a dominant microcrack as shown in Figure 15c. 



VII. Concluding Remarks 

The present methodology forms the basis for multiscale modeling of metallic fracture that is rooted at the 
atomistic level. Molecular dynamics simulations and a newly defined cohesive zone volume element were used to 
develop constitutive relationships describing interfacial deformation and fracture processes near a crack tip that were 
then recast within parameterized cohesive zone models and implemented as part of a finite element analysis. This 
methodology provides a physical, albeit limited to nanometer scales, constitutive relation of decohesion that can lead 
to the construction of a database of CZMs for studying the failure of polycrystal micro structures. In the present 
study the MD simulations considered nanoscale crack propagation along a common 299 grain boundary in 
aluminum and yielded significantly different traction-displacement relationships for the crack growing in ductile and 
brittle directions along the grain boundary. 

By adjusting the parameters of the bilinear CZM, a range of peak tractions and opening displacements were 
considered to represent variations in GB properties. These variations arise from the presence of vacancies, 
interstitials, and other defects in addition to the effect of the GB misorientation angle to the GB structure. The CZM 
then provides a numerically efficient means of considering a broad range of GBs incorporating mechanisms and 
configurations from the atomistic modeling that contribute to the variations in GB properties. 

Using the atomistically-derived traction-displacement results as a basis, two different parameterizations for a 
CZM were considered where the relative amount of elastic and plastic deformation was varied either through 
varying the fracture toughness, while keeping the peak traction constant, or through varying the peak traction, while 
maintaining a constant fracture toughness. In both cases, the increase of plasticity results in an increase in the 
process zone size at the crack tip. Preserving the fracture toughness resulted in a constant equilibrium crack length 
at a fixed strain, independent of the peak traction. The peak traction affects the initial strain at which the crack 
opens. The CZM was then used to illustrate the formation of a dominant crack in a simple aluminum polycrystal. 
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